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Photonuclear reactions using a laser Compton scattering (LCS) gamma source provide a new method for 
producing radioisotopes for medical applications. Compared with the conventional method, this method has the 
advantages of a high specific activity and less heat. Initiated by the Shanghai Laser Electron Gamma Source 
(SLEGS), we conducted a survey of potential photonuclear reactions, (y, n), (Y, p), and (vy, y’) whose cross- 
sections can be measured at SLEGS by summarizing the experimental progress. In general, the data are rare 
and occasionally inconsistent. Therefore, theoretical calculations are often used to evaluate the production of 
medical radioisotopes. Subsequently, we verified the model uncertainties of the widely used reaction code 
TALYS-1.96, using the experimental data of the '"Mo(y,n)”Mo, °Cu(y,n)“Cu, and ®Zn(y, p)*'Cu reac- 


tions. 


Keywords: Medical radioisotope, Photonuclear reaction, LCS, Cross section 


I. INTRODUCTION 


Radioisotopes are widely used for the diagnosis and ther- 
apy of different diseases owing to their nuclear-physical prop- 
erties [1, 2]. Diagnostic radioisotopes can provide functional 
and metabolic information for the early treatment of dis- 
eased regions that have not yet undergone morphological or 
structural changes. Currently, positron emission tomography 
(PET) [3] and single-photon emission computed tomography 
(SPECT) [4] are the two major diagnostic techniques. Short- 
lived ß *-emitting radioisotopes are often used for PET. Typi- 
cal examples include Ic, INI ; 15 O, and 18h which have half- 
lives of 20, 10, 2, and 110 min, respectively. For SPECT, 
y-ray emitting radioisotopes are frequently used, of which 
?"Tc with a half-life of 6 h [5] is the most common ra- 
dioisotope tracer. Therapeutic radioisotopes can be combined 
with targeted drugs to achieve the precise removal of small 
diseased regions without excessive doses to normal tissues. 
Therefore, radioisotopes that emit low-range highly ionizing 
radiation are of significant interest. The ß”-particle emit- 
ting radioisotopes (e.g. m 89Sr, Vy, oiy Lu, and 188Re), 
Auger electron cascades (e.g. Pd and T) and «-particle 
emitting radioisotopes (e.g. 21l At, p Ra, PAG and 
Th), with a high linear energy transfer in tissue, are suit- 
able for therapy. 

More than 40 million nuclear medicine procedures are per- 
formed annually, and the demand for radioisotopes is increas- 
ing by 5% annually [6]. Currently, the production of medi- 
cal radioisotopes relies primarily on nuclear reactors and cy- 
clotrons. Thermal neutron-induced fission produces neutron- 
rich radioisotopes in nuclear reactors. The advantage of this 
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method is the possibility of producing high levels of total 
and specific activities. However, the production of desired 
radioisotopes is accompanied by a considerable amount of 
long-lived radioactive waste, which raises numerous safety 
and security concerns. In addition, many nuclear reactors 
used for medical radioisotope production are more than 50 
years old and may be shut down in the near future [7]. For 
example, a vast majority of reactors producing ”Mo are ex- 
pected to shut down by 2030. Compared with reactors, cy- 
clotrons typically produce neutron-deficient radioisotopes by 
charged particle reactions accompanied by less radioactive 
waste. In addition, they have the advantage of being more 
compact, allowing their placement near hospitals, and mak- 
ing them a useful tool for producing radioisotopes with short 
half-lives that range from several minutes to hours. 

Another promising method for radioisotope production 
involves photonuclear reactions, mainly including (Y,n), 
(y, p), and (y,y’). This is considered an alternative to ra- 
dioisotope production in reactors and cyclotrons or the only 
method of production for certain radioisotopes. (y,7) reac- 
tions can produce ß*-emitting radioisotopes for PET imag- 
ing, such as Ic, DN, 150, and "®F. These radioisotopes are 
produced by cyclotron-based (p, n) or (p, y) reactions. (y, p) 
reactions are suitable for producing f -emitting radioiso- 
topes that are currently mostly produced in reactors. High- 
intensity y beams are required to obtain adequate yield using 
this method. With the development of laser Compton scat- 
tering (LCS) gamma sources, the production of radioisotopes 
via photonuclear reactions has attracted considerable atten- 
tion [8-17]. 

In contrast to conventional bremsstrahlung gamma sources 
based on electron linear accelerators, LCS gamma sources 
have the advantages of a high photon flux and excellent 
monochromaticity. The advantages of photonuclear reactions 
using LCS gamma sources for radioisotope production are as 
follows: 


* Differing from bremsstrahlung gamma sources, the 
LCS gamma sources significantly reduce the heat per 
useful reaction rate as its y-rays in the energy range of 


interest are not accompanied by an intense low-energy 
tail. Moreover, the LCS gamma sources have the abil- 
ity to selectively tune photons to energies of interest. 
A high specific activity can be achieved by matching 
the photon energy to the giant dipole resonance (GDR) 
peak. 


Differing from the charged-particle induced reactions, 
the photonuclear reactions generate less heat as their 
energy deposition in targets via gamma-matter inter- 
actions is significantly smaller. Thus, the target can 
be thicker and requires considerably less cooling [9]. 
Moreover, the photonuclear reactions have the capabil- 
ity to simultaneously irradiate multiple targets. In this 
manner, it can maximize the utilization of y beams and 
produce multiple radioisotopes simultaneously [8]. 


The Shanghai Laser Electron Gamma Source (SLEGS) 
employs the LCS technique. It can generate y beams 
ranging from 0.25-21.7 MeV with a full-spectrum flux of 
10°-107 s-! [18] and the best possible bandwidth of y beams 
of 5%-15% after passing through a dual collimation sys- 
tem [19]. Using collimation technology, the size of an y beam 
can be continuously adjusted to within $25 mm [20]. The 
monochromaticity and high intensity of y beams combined 
with detector spectrometers can be used to measure the pho- 
tonuclear reaction cross-section [21—24]. One of the main 
topics of SLEGS is the measurement of the key photonuclear 
reaction cross-section relevant for medical radioisotope inves- 
tigations, providing crucial nuclear data for the photonuclear 
method. 

In this study, we aim to estimate the production of med- 
ical radioisotopes in the energy range of 0.25-21.7 MeV in 
the SLEGS domain. The remainder of this paper is organized 
as follows. In Sect. II, we summarize the radioisotopes that 
can be produced via (y, n), (y, p), and (y, y') reactions and 
discuss their experimental feasibility. Owing to the scarcity 
of cross-sectional data regarding photonuclear reactions, the- 
oretical calculations are often used to assess the yields of 
medical radioisotopes. In Sect. III, the cross-sectional data 
of the !°Mo(y,n), 9Cu(y, n), and Zn(y, p) reactions in 
the SLEGS energy region were investigated to produce Mo, 
Cu, and "Cu radioisotopes. Finally, we provide a summary 
in Sect. IV. 


Il. MEDICAL RADIOISOTOPES PRODUCED IN 
PHOTONUCLEAR REACTIONS 


In this section, we discuss potential (y,n), (y,p), and 
(y, y’) reactions for the production of medical radioisotopes. 
Crucial considerations include the natural abundance of the 
target isotopes, half-life of the radioisotopes, and available 
experimental data. Generally, it is preferable to have a target 
isotope with a dominant abundance (e.g., , greater than 10%), 
and a suitable half-life (several tens of minutes to days) for 
the produced radioisotope. However, the lack of experimen- 
tal data may limit the medical applications of radioisotopes. 


Different methods are available for measuring photonu- 
clear reactions, including in-beam and offline measurements. 
Four 7 ?He neutron detection arrays, charged particle detec- 
tors, and nuclear resonance fluorescence spectrometers are 
commonly used for in-beam measurements of (y, n), (y, p), 
and (y, y^) reactions. Offline measurement is suitable for de- 
termining all three types of photonuclear reactions only in 
cases where the residual decay has an unambiguous Y-ray 
and the lifetime of the residual nucleus is within a range of 
minutes to hours. 


A. (y,n) reactions 


Radioisotopes can be most efficiently produced by excit- 
ing the target nuclide into a GDR using photons. The GDR is 
characterized by a large peak cross-section and broad width, 
resulting in a large integral cross-section. The shape of the 
GDR follows a Lorentzian distribution with a spreading width 
of approximately 5 MeV [25-27]. Generally, a consider- 
able portion of the total photonuclear cross-section is known 
to originate from the (y,n) reactions, which increases as 
the charge number of the target nuclide increases, making 
the following chemical separation easier. For example, the 
MMo(y, n) reaction accounts for more than 70% of the to- 
tal photonuclear reactions within the energy range of 8.4-15 
MeV. Therefore, medium- and heavy-mass nuclei are ideal 
for production via (y, n) reactions [28]. 

Table 1 summarizes the 15 medical radioisotopes produced 
by (y,n) reactions [8, 14, 29-31]. Excluding Mo, the 
natural abundance of these reaction targets exceeds 10%. A 
higher natural abundance results in a higher yield and lower 
impurity content. The seventh column indicates the current 
status of the EXFOR database [32]. The experimental data 
for only nine radioisotopes are presented in the Appendix. 
Among these, the PR(y, n) reaction has only one data set 
that was measured in the 1960s. The data for the °°Zr(y, n), 
00Mo(y,n), and '"Re(y, n) reactions do not cover the en- 
tire GDR energy region. There is a significant discrepancy 
between the two sets of data for the Scu(y, n) reaction: It is 
possible to determine the relevant cross-sections at SLEGS. 


B. (y,p) reactions 


Despite the nucleus being excited beyond the proton sep- 
aration energy of the photons, it does not necessarily lose a 
proton. Only for excitations well beyond the proton separa- 
tion energy, the proton can acquire sufficient kinetic energy 
to effectively cross the Coulomb barrier. However, in these 
cases, the excitation energy typically exceeds the separation 
energies of one or two neutrons. Therefore, the neutron emis- 
sion channel competes with the proton emission channel. For 
light and certain medium-mass nuclei, the cross-sections of 
the (y, p) reactions are comparable to and occasionally ex- 
ceed those of the (y, n) reactions, owing to the shell structure 
of the nuclei. Certain heavier D^ emitters for radionuclide 
therapy can also be produced by (y, p) reactions, but the in- 


TABLE 1. Medical radioisotope ^X production via the ^+! X (y, n)^ X reactions [8, 14, 29-31]. The natural abundance of ^*! X is 
provided. The produced medical radioisotopes ^ X are characterized by the half-lives T; /2 and decay modes. Here, %e represents the proba- 
bilities of nuclear decay via electron capture or ß* decay. %ß ^ represents the probabilities of ^ decay. E; represents the threshold energy 
of the reaction. The seventh column indicates whether the ^*! X (y, n)^X reactions have experimental data in the EXFOR database [32]. a 
and b indicate that the radioisotope can be studied by in-beam and off-line measurements, respectively. The last column presents the probable 
y-rays of interest in the offline measurement. The data are obtained from Ref. [33]. 


AFTY Nat. abu. (%) ay T2 Decay mode Fin (MeV) EXFOR data Medical application Exp. method y-ray (keV) 
Pc 98.93 lC 20364m 96e = 100 18.7 Yes PET a;b 511 
N 99.63 PN 9965m 96e = 100 10.5 Yes PET a;b 511 
169 99.757 PO 12224s %e=100 15.6 Yes PET a z 
IE 100 IF 109.77m 96e — 100 10.4 Yes PET a;b 511 
BSc 100 “sc  397h 96e = 100 11.3 No PET a;b 511;1157 
8Cu 69.15 Cu 9673m %e=100 10.8 Yes PET a;b 511;1157 
Scu 30.85 “Cu  127h %e=61.5 9.91 Yes PET;Radiotherapy a;b 511 

%ß = 38.5 
$9 Ga 60.108 Ga 6771m 96e = 100 10.3 No PET a;b 511;1077 
907t 51.45 Zr 78.41h 96e = 100 11.9 Yes PET a - 

10Mo 9.82 Mo 65976h %ßT = 100 8.29 Yes SPECT a;b 740;778 
104pq 11.14 pq 17.0d %e = 100 10.0 No Radiotherapy a s 
\eORr 33.503 Rr 1036h 96e — 100 8.47 No Radiotherapy a - 
UOge 14.91 Fr 94d %B~ = 100 7.25 No Radiotherapy a : 

WRe 62.60 "Re 37d %B~ = 92.53 7.36 Yes SPECT;Radiotherapy a - 
He = 7.47 

Tr 62.70 Tr 738d APß = 95.24 7.77 No Radiotherapy a : 
96e — 4.76 


creasing Coulomb barrier leads to the production of small 
cross-sections. (y,p) reactions result in the daughter and 
parent isotopes being chemically different. Several advanced 
chemical-separation techniques have been developed for this 
purpose. For example, a simple and reproducible three-ion 
exchange matrix approach was developed to separate Cu 
from Zn [34]. Separation techniques exist for other pairs, 
such as Ti/Sc, Zr/Y, and Hf/Lu [35-37]. 

The potential (y, p) reactions used for nuclear medicine are 
listed in Table 2. However, the experimental data regarding 
this topic are scarce. Only limited data for ^K, °’Cu, and 
"Tu has been obtained by using the bremsstrahlung gamma 
source, as detailed in the Appendix. No error bars are shown 
for the ?K data, whereas relatively large errors ranging from 
typically 10% to 100% for "Cu and '’’Lu have been demon- 
strated. 


C. (y,y) reactions 


A nuclear isomer is the metastable state of an atomic nu- 
cleus, in which one or more nucleons (protons or neutrons) 
occupy higher energy levels than the ground state of the same 
nucleus. The lifetime and excitation energy are two important 
properties of nuclear isomers. The excitation energy between 
the isomer and ground state is characteristic and can be used 
to identify nuclides. Nuclear isomers can be deexcited by 
the emission of y-rays and/or the conversion of electrons to 
the ground state. The y-emitting isomers can be used as ra- 
dioactive labels for SPECT imaging. Nuclear isomers emit- 
ting low-energy Auger electrons are potential radioisotopes 


for targeted therapy. 

Conventional production methods, such as (n, y) reactions, 
have relatively low yields because the dominant part of the 
production proceeds directly to the nuclear ground state with 
a spin closer to that of the target isotope. Using monochro- 
matic and small-bandwidth LCS photons, transitions from a 
stable or long-lived nuclear ground state to higher energy lev- 
els can be selectively excited. Such levels serve as gateway 
states, which then partially decay to the isomeric state directly 
or via a cascade. The (y, y’) reaction is equivalent to storing 
the energy of an incident photon in an isotope that acts as a 
container. 

Table 3 lists the six isomers that can be produced by (vy, y) 
reactions for medical applications. Excluding !'""Sn and 
176m u, the experimental data for the other radioisotopes are 
detailed in the Appendix. The dataset of "Rp did not have 
error bars, and two different results were reported from 6 to 
23 MeV. For !"In, two cross-sectional datasets differ by 
more than 300 times at 8 MeV. For ''””In, there were sig- 
nificant discrepancies between multiple datasets. The data of 
P5"Pt have a large relative error as high as 100% at 3.5, 7.5, 
and 8 MeV. 


IH. EVALUATION OF CROSS-SECTIONS FOR ?MO, 
“CU, AND "CU PRODUCTION 


The cross-section data were used to calculate the expected 
yield of the radioisotopes for a given thickness and enrich- 
ment of the target material, and to determine the optimum 
energy range for the production of the desired radioisotope 


TABLE 2. Same as Table 1 but for medical radioisotopes AX production via the Atly(y, p)^X reactions [8, 38—40]. Nat. abu. indicate the 
natural abundance of ^*^ Y. The seventh column indicates whether experimental data regarding the AHy(y, p)^X reactions is available in 
the EXFOR database [32]. 


+1y Nat. abu. (96) X Ty/2 Decay mode E;n (MeV) EXFOR data Medical application Exp. method y-ray (keV) 
Ca 2.09 “K 223h %ß” — 100 12.1 Yes Radiotherapy a;b 373;617 
‘STi 73.72 "$c  335d %ßT = 100 11.4 No Radiotherapy a : 
Ni 68.077 "Co 27174d %e = 100 8.17 No Radiotherapy a S 
Zn 18.45 Cu 61.83h %ßT = 100 9.97 Yes Radiotherapy a;b 861;2502 
u^ 11.22 VY  64053h 968- = 100 8.68 No Radiotherapy a : 
!06pd 27.33 Rh 3536h %ßBT = 100 9.34 No Radiotherapy a;b 319 
IPS 0.97 Hiin 2805d %e= 100 7.55 No SPECT;Radiotherapy a - 
132Xe 26.9 BII 8.025d %ß” = 100 9.12 No Radiotherapy a - 
Dy 25.475 "Tb 689d %ßBT = 100 8.00 No SPECT;Radiotherapy a - 
Er 22.869 ‘Ho 26.824h 968- = 100 7.50 No SPECT;Radiotherapy a : 
ef 27.28 Lu 6647d %ß = 100 7.34 Yes SPECT;Radiotherapy a;b 208 
1890s 16.15 Re 17h —968- — 100 725 No SPECT;Radiotherapy a : 


TABLE 3. Same as Table | but for medical isomers ^"" X production via the ^ X (v. Y! yim reactions [8, 13]. The isomeric states 4” X 
are characterized by the energy Em, half-lives 1/2: and decay modes. Here, %IT represents the probabilities of nuclear decay via isomeric 


transitions. Ey, equals Em. The seventh column indicates whether the ^ X (y, y)?" X reactions have experimental data in the EXFOR 
database [32]. ^ X (y, y)?" x reactions can be investigated by off-line measurements where the y-ray energy of interest is equal to Em. 


X Nat. abu. (46) X Ti /2 Em/Eın (keV) Decay mode EXFOR data Medical application 
Rh 100 de 


Rh 56.114 m 39.8 %IT = 100 Yes Radiotherapy 
"In 4.29 "By 99.476 m 391.7 %IT = 100 Yes SPECT;Radiotherapy 
1I 95.71 Inm 4.486 h 3362 %IT — 95 Yes SPECT;Radiotherapy 
AB =5 
Sn 7.68 men 13.76 d 314.6 %IT = 100 No Radiotherapy 
Lu 2.599 i E 3.664 h 122.8 %ß = 99.90 No Radiotherapy 
Ae = 0.09 
I95pr 33.78 I93mpt 4.01 d 259.3 %IT = 100 Yes Radiotherapy 


and the level of radioisotopic impurities. However, the slow 
development of gamma sources and small cross-sections of 
photonuclear reactions resulted in a relatively small num- 
ber of experimental studies regarding this cross-sectional 
data. There is often scattering between different experimental 
datasets. As a result, the cross-sections predicted by the the- 
oretical models play a key role in the production of medical 
radioisotopes. 


The TALYS code [41] offers a unified approach for calcu- 
lating nuclear reactions that involve neutrons, photons, pro- 
tons, deuterons, tritons, He, and « particles in the keV-200 
MeV energy range and for target nuclides with masses of 5 
and 339. The code outputs a set of reaction data, for example 
cross-sections, energy spectra, and angular distributions of 
the emitted particles, etc.. Numerous studies have tested the 
TALYS code and demonstrated that it has reliable predictive 
power in terms of the nuclear reaction calculations [42-44]. 
In the TALYS code, the decay of the compound nuclear state 
from the photonuclear reaction was treated using the Hauser- 
Fenshbach (HF) statistical model [45]. The nuclear-level den- 
sity and y strength function were the main input parameters. 
In this study, we adapted the TALYS-1.96 code. It imple- 
mented six different nuclear level density (NLD) models and 
nine different models of y strength functions (ySF), as shown 
in Table 4. The NLD and ySF models employ phenomeno- 


logical and microscopic approaches, respectively. 

In this section, we select the Mo, Cu, and Cu ra- 
dioisotopes as examples to compare the model calculations 
with the experimental data. ””Tc decayed by ”Mo is the 
most frequently used radioisotope in nuclear medicine. Zn 
is of interest owing to its applications in the production 
of the medical radioisotope Cu; 'Cu and “Cu can form 
a "matched pair" [46, 47]. Therapeutic "Cu, along with 
positron-emitting Cu, can measure the uptake kinetics in 
an organ of a patient by PET imaging, allowing for a precise 
dosimetric calculation. We used the default values of NLD 
and the ySF models in the following calculations. 


A. 9 Mo/?"Tc 


With half-lives of 6 h and 140 keV y-rays, ??"Tc is 
nearly ideal for SPECT imaging. The most common method 
to obtain ?""Tc is by elution from ??Mo generators. The 
10Mo(y,n) reaction has a sufficient potential to produce 
Mo. The excitation functions are presented in Fig. 1, which 
covers the entire GDR energy range. The cross-sections 
of 'Mo(y, n)?Mo were experimentally measured by Ut- 
sunomiya et al. [48], Crasta et al. [49], and Ejiri et al. [14]. 
The gamma sources used in the experiments were LCS and 


TABLE 4. The descriptions of NLD and SF models in the TALYS-1.96 code 


Input parameter 
Idmodel 1 (default) 
Idmodel 2 
Idmodel 3 
Idmodel 4 
Idmodel 5 
Idmodel 6 
strength 1 (default for incident neutrons) 
strength 2 (default for other incident particles) 
strength 3 
strength 4 
strength 5 
strength 6 
strength 7 
strength 8 
strength 9 


bremsstrahlung. All the experimental results were distributed 
in the low-energy region of the GDR and did not exceed the 
peak. Fig. 1(a) demonstrates that the shape and value of the 
excitation function are highly sensitive to the choice of the 
ySF model, particularly at the peak position of GDR. Ex- 
cluding Strength 2 and 3, the remaining YSF models yield 
relatively similar results that are close to the experimental re- 
sults. Different ySF models produce similar predictions for 
energies above the GDR peak. Overall, the cross-sections es- 
timated using Strength 8 achieve the best agreement with the 
experimental data. As shown in Fig. 1(b), the calculation re- 
sults of all the NLD models below the peak are consistent but 
significantly higher than the experimental data, whereas those 
above the peak are significantly different. The data predicted 
by the models at energies above 14 MeV are highly valuable. 
New experiments can be conducted based on SLEGS to verify 
the model predictions. 


B. “Cu 


As listed in Table 1, the decay characteristics of Cu ren- 
der it useful for nuclear medicine [50]. It combines PET di- 
agnostic capabilities with those of radiotherapy, with average 
electron emissions of 190 keV. Currently, “Cu is mainly pro- 
duced by small cyclotrons through Ni (p, n) reactions [51— 
54]. Alternative production using 6 Cu(y, n) does not require 
rare or expensive “Ni targets and simplifies the chemical sep- 
aration step. The experimental data for the cross-sections 
of the SCu(y, n) Cu reaction are presented in Fig. 2 along 
with the TALYS. The measurements by Katz et al. [55] and 
Antonov et al. [56] were performed using bremsstrahlung 
gamma sources. However, they differed significantly in terms 
of their shape and value. Coote et al. [57] obtained a sin- 
gle cross-section at 17.6 MeV using monochromatic y-rays 
from the "Li(p, y)?Be reaction. The cross-section is consis- 
tent with the calculation but lower than other experimental 
data by 80 to 120 mb. The cross-sections of 9Cu(y, n) “Cu 
can significantly vary with the choice of the ySF models but 
not the NLD models. However, the discrepancies between the 


Detailed description 
Constant Temperature + Fermi gas model 
Back-shifted Fermi gas model 
Generalised Superfluid model 


Skyrme-Hartree-Fock-Bogoluybov level densities from numerical tables 
Gogny-Hartree-Fock-Bogoluybov (GHFB) level densities from numerical tables 
Temperature-dependent GHFB level densities from numerical tables 


Kopecky-Uhl generalized Lorentzian 
Brink-Axel Lorentzian 


Hartree-Fock + Bardeen-Cooper-Schrieffer (BCS) tables 


Hartree-Fock-Bogoliubov (HFB) tables 
Goriely’s hybrid model 
Goriely Temperature-dependent HFB model 


Temperature-dependent Relativistic Mean Field (RMF) model 
Gogny D1M HFB + Quasiparticle Random Phase Approximation (QRPA) 


Simplified Modified Lorentzian (SMLO) model 
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Fig. 1. (Color online) "Mo(y, n)”Mo reaction cross sections. 


experiments and models cannot be compensated for by vary- 
ing the y SF models. Future experiments regarding the Cu 
photoneutron reaction using a monochromatic LCS gamma 
source will be useful for resolving these discrepancies and 
provide assurance for its medical applications. 


EA LELEI a a vn a a u vv a 994 99 0909 00 09 
200 } Bs E 
——— strength] o 
180 - — strength? E | 
——— strength3 
—— strength4 H o 
160 F —— strength5 |] 
—— strength6 A A 
© 140 L —— strength7 A H Bn -| 
g€ — strength8 hn 
= strength9 o 
g Dor © 1990-Antonov 4 B j 
.g -e-! 1961-Coote 
Ə 100+ 4 1951-Katz 
[5] 
n 
2 807 
g 
O 60- 
40 r 
20 F 
0 E 
| | | J | | J | | | | | 
T T T T T T T T T T T T 
200 F Be al 
o 
180 L — Idmodell = J 
—— |dmodel2 n 
—— — ]dmodel3 [m] 
160 F — Idmodel4 J 
= — — Idmodel5 n^ A 
2 140+ —— Idmodel6 ya B 4 
g Li 1990-Antonov N 
— He 1961-Coote o n 
g 1207 4 1951-Katz 2 7 
[s 
B 
Q 
[5] 
a 
n 
a 
fe) 
= 
U 
l | l | J i l | l | l | 


10 1] 12 13 14 15 16 17 18 19 20 21 
Energy [MeV] 


Fig. 2. (Color online) 9Cu(y, n)™Cu reaction cross sections. 


C. "Cu 


6’Cu is a promising f -particle emitter with an average en- 
ergy of 141 keV for the targeted radiotherapies. The size of 
these particles in the tissue was of the same order as the cell 
diameter. This reduces the unwanted dose burden on patients 
[58]. With a half-life of 61.83 h and low energy y emissions 
(91.266 keV, 7%; 93.311 keV, 16.1%; 184.577 keV, 48.7%), 
6’Cu can provide a long therapeutic effect. "Cu and its stable 
daughter °’Zn are nontoxic to the body, and both Cu and Zn 
are prevalent trace elements in the body. Along with the PET 
imaging radioisotope Cu, it forms a “matched pair.” How- 
ever, the widespread clinical use of °’Cu-based radiopharma- 
ceuticals has been limited by their availability, quantity, and 
quality. Experiments [59-63] have shown that the °’Zn(y, p) 
reaction has the potential to produce sufficient quantities of 
6’Cu with a sufficient purity for medical use. Moreover, var- 
ious Cu and Zn separation methods have been employed in 
radiochemical processing [34, 64-66]. 

Figure 3 presents the excitation function of the 
87n(y, p)Cu reaction. Experimental data regarding 
this reaction are scarce, with only one group of two points in 
the energy region of SLEGS [67]. The data are significantly 
inconsistent with any calculation and hardly constrain NLD 


and ySF models. In contrast to the (y, n) reaction of several 
hundred mb, the (y, p) reaction is typically several mb in 
the medium-heavy nuclei. Additionally, protons can easily 
stop at the target and are difficult to detect. If the product 
nuclei of the (y, p) reaction are unstable and have a moderate 
half-life, the cross-sections can be determined by offline 
measurements of the characteristic y-rays during decay. For 
the 687 n(y, p) reaction, an offline measurement technique 
can be employed to determine its cross-sections. 
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Fig. 3. (Color online) ®Zn(y, p) S7Cu reaction cross sections. 


IV. SUMMARY 


In this study, we summarize the medical radioisotopes that 
can be produced by (y, n), (y, p), and (y, y’) reactions, along 
with the experimental data for these production pathways. We 
investigated the photonuclear reaction cross-sections for the 
production of "Mo, “Cu, and Cu. The experimental data 
for Mo and °’Cu did not cover the entire GDR energy re- 
gion, whereas the data for Cu exhibited a significant dis- 
crepancy. To better constrain the NLD and ySF models in 
TALYS, additional experimental measurements at SLEGS in 
the energy region of GDR are necessary in the future. 


To compare the productivity of the photonuclear method 
with that of traditional methods, we consider the produc- 
tion of ”Mo/””Tc as an example. The highest flux cur- 
rently available for SLEGS is 10" s^!. We adopted the cal- 
culated results by the TALYS-1.96 code using the Strength 
8 model for the yield calculations, which provided a max- 
imum of 145 mb of the '°°Mo(y, n) reaction cross-section 
at the incident photon energy of 14.4 MeV. After irradi- 
ating a '°°MoO; target with a radius of 2 mm and a thick- 
ness of | cm for 5 times the half-life of Mo, the saturation 
specific activity of ”Mo and ””Te can reach 4.76 and 4.74 
uCi/g, respectively. The use of stacked targets can further en- 
hance production. When the energy regions corresponding to 
the maximum cross-sections of the targeted nuclear reactions 
are close, multiple radioisotopes can be simultaneously pro- 
duced. For example, the Cu, “Cu, and ®°Zr radioisotopes 
can be produced simultaneously when the incident photon en- 
ergy is approximately 17 MeV. 


Currently, Mo is primarily produced by the (n, f) re- 
action in a high-flux reactor using an enriched "U target. 
The specific activity of Mo can reach 185 TBq/g (5000 
Ci/g) [6]. However, most reactors will gradually shut down 
by 2030 [68]. Medical cyclotrons can be used to directly 
generate "Tc via the '°°Mo(p, 2n) reaction. According to 
the experimental data, an enriched Mo target irradiated 
with a 16.5 MeV proton beam at 130 uA for 6 h yielded 
116 GBq/g (3.13 Ci/g) of "Tc [69]. Recently, a subcrit- 
ical Mo production system was developed, which is driven 
by an accelerator-based deuterium-deuterium (D-D) neutron 
source. The D-D fusion reaction generates neutrons that ir- 
radiate a low-enriched uranium solution and induce fission in 
2350, The system can generate 47.8 mCi /g Mo for a stable 
24 h operation with a neutron intensity of 1x 10!4 n/s [70]. 
Bremsstrahlung gamma sources based on linear electron ac- 
celerators are often used to produce ??Mo radioisotopes. Us- 
inga INO, target irradiate with a 35 MeV electron beam 
at 100 uA for 20 h; the specific activity can achieve 4.4 
GBq/g (119 mCi/g) [71]. NorthStar Medical Radioisotopes 
developed a "Mo production system using an electron lin- 
ear accelerator. Using this system based on (y, n) reactions, 
approximately 3096 more ??Mo is produced per gram of the 
target material compared with the traditional neutron capture 
route [72, 73]. 


In comparison, the specific activity produced by the pho- 
tonuclear method based on SLEGS was low. With a fur- 
ther increase in the flux of y beams, from 10" s^! to 10% 
s-! [20], the specific activity of Mo can reach up to 4x 10? 
Ci/g. For the 2-day protocol, the activity of the ””Tc- 
labelled tracers required for one myocardial perfusion imag- 
ing was 24 mCi [74]. The total yield of stacked targets in 
one year can provide 400,000 myocardial perfusion images. 
As the flux of SLEGS increases in the future, this method is 
promising for producing medical radioisotopes and will be- 
come feasible in China. 


APPENDIX 


Experimental data for the production of medical radioiso- 
topes by (y,n), (y,p), and (y, y') reactions are summa- 
rized. Table 5 lists relevant information on these reactions. 
All the data and information were obtained from the EXFOR 
database [32]. 


TABLE 5. Information regarding experimental data for the production of medical radioisotopes by photonuclear reactions 


Target Reaction Product Year Gamma source Reference 


od (y, n) Hc 1951 Bremsstrahlung [55] 
1957 Bremsstrahlung [75] 

1959 ?*H(p, y)*He [76] 

1961 Bremsstrahlung [77] 

1962 ?H(p, y)^He [78] 

1966 Bremsstrahlung [79] 

1966 Bremsstrahlung [80] 

1966 ?*H(p, y)*He [81] 

"N (y,n) BN 1960 Bremsstrahlung [82] 
1987 Bremsstrahlung [83] 

160* (y,n) 150 1962 ?H(p, y)^He [78] 
1966 Bremsstrahlung [84] 

1970 Bremsstrahlung [85] 

1985 Bremsstrahlung [86] 

1991 Bremsstrahlung [87] 

PE (y, n) Ip 1962 ?H(p, y)^He [78] 
cu (y, n) Cu 1951 Bremsstrahlung [88] 
1951 Bremsstrahlung [55] 

1954 Bremsstrahlung [89] 

1955 Bremsstrahlung [90] 

1959 "Li(p, y) *Be [91] 

1960 "Li(p, y) *Be [92] 

1961 "Li(p, y)*Be [57] 

1962 ?H(p, y)^He [78] 

1968 In-flight positron annihilation [93] 

1968 Bremsstrahlung [94] 

1979 In-flight positron annihilation [95] 

1984 "Li(p, y)*Be [96] 

1984 Bremsstrahlung [97] 

2012 Bremsstrahlung [98] 

Cu (vn) cu 1951 Bremsstrahlung [55] 
1961 "Li(p, y)’Be [57] 

1990 Bremsstrahlung [56] 

Zr (y, n) 97r 1978 Bremsstrahlung [99] 
2010 Bremsstrahlung [100] 

2019 LCS [101] 

10Mo (y, n) Mo 1978 Bremsstrahlung [99] 
2011 Bremsstrahlung [49] 

2011 LCS [14] 

2013 LCS [48] 

187Re (y,n) Re 2005 LCS [102] 
“Ca (y, p) BK 1978 Bremsstrahlung [99] 
Zn (y, p) Cu 1957 Bremsstrahlung [67] 
2021 Bremsstrahlung [56] 

Ug (v, p) "py 2020 Bremsstrahlung [103] 
IRh Cv, Y!) Impp 1961 Bremsstrahlung [104] 
Mig (y, Y) Ben 2006 Bremsstrahlung [105] 
2017 Bremsstrahlung [106] 

5m Cy. Y) Hn 1957 Bremsstrahlung [107] 
1993 Bremsstrahlung [108] 

2001 Bremsstrahlung [109] 

2006 Bremsstrahlung [110] 

2017 Bremsstrahlung [106] 

2018 LCS [111] 

Spr (v, y’) I95mpr 2006 Bremsstrahlung [105] 


* Experimental data within a specific energy range, instead of all of experimental data 
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Fig. 6. (Color online) 'éO(y, n.) ^O reaction cross sections. Fig. 9. 8Cu(y, n)“Cu reaction cross sections. 
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Fig. 10. °Zr(-y,n)*Zr reaction cross sections. 
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Fig. 11. '°Mo(y,n)”Mo reaction cross sections. 
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